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Isospin effects of the critical phenomena were studied via Xe isotopes in the frame of lattice 
gas model. All the critical temperatures for four Xe isotopes are close to 5.5 MeV at the same 
freeze-out density of about 0.39 po. The critical values of power law parameter of mass distribution, 
mean multiplicity of intermediate mass fragments (IMF), information entropy and Campi's second 
moment show minor dependence on the isospin at the critical point. 



O i With the development of the accelerator and radioactive beam technique, a lot of new experiments using the 
radioactive beams with large neutron or proton excess become possible. The degree of freedom of isospin of nuclear 
matter is becoming important for research. It offers the possibility to study the properties of nuclear matter in 
the range from symmetrical nuclear matter to pure neutron matter. Some theoretical investigations to the equation 
of state, chemical and mechanical instabilites as well as liquid-gas phase transition for isospin asymmetrical nuclear 
matter were performed already. In addition, the isospin dependent nucleon - nucleon cross section is also an important 
subject due to its significant effects on the dynamical process of heavy ion reactions induced by radioactive beams. 
^sO , Some new phenomena stemmed from the isospin have been revealed. For examples, the isospin dependences of the 
1 preequilibrium nucleon emission, the nuclear stopping, the nuclear collective flow, total reaction cross section, radii 
of neutron-rich nuclei and subthreshold pion production have been studied by several groups However, more 

experimental and theoretical studies are still needed for understanding the isospin physics. As a trial, the isospin 
' effects were investigated with the lattice gas model in this letter. 

, The lattice gas model of Lee and Yang |I| , in which the grancanonical partition function of a gas with one type 
of atoms is mapped into the canonical ensemble of an Ising model for spin 1/2 particles, has uccessfully described 
the liquid-gas phase transition for atomic system. The same model has already been applied to nuclear physics for 
1 1 ■ isospin symmetrical systems in the grancanonical ensemble Q with an approximate sampling fl o f the canonical 
ensemble Jz| — P~0|| , and also for isospin asymmetrical nuclear matter in the mean field approximation JLl| . In this model, 
A nucleons with an occupation number s which is defined as s = 1 (-1) for a proton (neutron) or s = for a vacancy, 
are placed in the L sites of lattice. Nucleons in the nearest neighbouring sites have interaction with an energy e SiS . 
The hamiltonian is written by 
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The interaction constant e SiSi is fixed to reproduce the binding energy of the nuclei, e nn ^ P p = e-i-i.n = 0. MeV, 
£pn,np = £1-1,-11 = -5.33 MeV. We use a three-dimension cubic lattice L with a size 1, a number of nucleons A = N + Z 
and a temperature T. The freeze-out density of disassembling system is pj — j^po where po is the normal nucleon 
density. The disassembly of the system is to be calulated at pf, beyond which nucleons are too far apart to interact. 
A nucleons are put in L cubes by Monte Carlo sampling using the Metropolis algorithm [H . Once the nucleons have 
been placed, their momentum is generated by a Monte Carlo sampling of Maxwell Boltzmann distribution. Various 
observables can be calculated in a straightforward fashion. 

One of the basic measurable quantities is the distribution of fragment mass. In this lattice gas model, two neigh- 
boring nucleons are viewed to be in the same fragment if their relative kinetic energy is insufficient to overcome the 
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attractive bond: P 2 /2/x + e SiSj < 0. This method is similiar to the so-called Coniglio-Klein's prescription jlj|. In this 
letter, we use the above condition to construct the fragments and their distributions. 

We chose several isotopes of Xe as examples of the study of isospin effects in the lattice gas model. Their isospin 
parameter {^-^-) is 0.11, 0.16, 0.21 and 0.26 for 122 Xe, 129 Xe, 137 Xe and 146 Xe, respectively. The freeze-out density 
p f has been chosen to be close to 0.39 po, extracted from the analysis of Ar + Sc |tJ and 35 C1 + Au and 70 Ge + Ti 
[14 1 with the same model. There is also good support from experiment that the value of pf is significantly below 0.5po 
[ 1 5 1 . We use the 343 cubic lattice with size of 7 which results that the freeze-out density Pf/po of 122 > 129 > 137 446x e 
is 0.36, 0.38, 0.40, and 0.43, respectively. The other input parameter is the temperature, we perform the calculation 
from 4 to 7 MeV. For each isotope 1000 events are accumulated at each temperature. 

Fig.l shows the mass distribution of fragments at T = 4, 5, 6 and 7 MeV for 129 Xe. Clearly the disassembling 
mechanism evolves with the nuclear temperature. A few light particles and fragments are emitted and the big 
residue reserves at T = 4 MeV which indicates typical evaporation mechanism. With the increasing temperature, 
the shoulder of mass distribution occurs due to the competition between the fragmentation and the evaporation. 
This shoulder disappears and the mass distribution becomes power law shape at T = 6 MeV, corresponding to the 
multifragmcntation region. When the temperature becomes much higher, the mass distribution becomes steeper 
indicating that the disassembling process becomes more violent. The power law fit, Y(A) oc A~ T , for these mass 
distribution can be introduced here. It has already been observed that a minimum of power law parameter T m i n exists 
for most systems if the critical behavior takes place. The lines in Fig.l represent the power law fit. Fig. 2 displays 
the several physical quantities as a function of temperature for Xe nuclei with the different isospin. The minimums 
of t parameters in Fig. 2a locate closely at 5.5 MeV for all the systems, which illustrates its minor dependence on 
the isospin. In other words, there is a universal mass distribution regardless of the size of disassembling source when 
the critical phenomenon takes place. However, the r parameters show different values outside the critical region for 
nuclei with different isospin, eg., r decreases with isospin when T > 5.5 MeV (multifragmentation region). Similiarly, 
the mean multiplicity of intermediate mass fragment Njmf-> defined as the number of fragments with 3< Z <16 here, 
has analogous characters in Fig. 2b |l^,|l8|. There are the maximums for Xe systems near to 5.5 MeV. When the 
temperature becomes higher, the larger the source, the higher the Njmf- 

Fig. 2c plots the information entropy H as a function of temperature for Xe isotopes. The information entropy was 
introduced by Shannon in information theory first M] . It is defined as 

H = -Y^Pilnfa), (2) 

where pi is the probability having "i" produced particles in each event, the sum is taken over all multiplicities of 
products from the disassembling system. H reflects the capacity of the information or the extent of disorder. We 
introduce this entropy into the nuclear disassembly here. As expected, the entropy H reveals the peak close to 5.5 
MeV for all isotopes. These peaks indicate that the opening of the phase space and the number of the states at the 
critical point is the largest. In the other words, the systems at the critical point have the largest fluctuation which 
leads to the largest disorder. After the critical point, the entropy H increases with the isospin and/or the source size. 

In Fig. 2d we give the temperature dependences of Campi's second moment of the mass distribution [fL6|, which is 
defined as 

b 2 7 , (3) 

where is the number of clusters with Ai nucleons and the sum excludes the largest cluster A max , A is the 

mass of the system. At the percolation point S2 diverges in an infinite system and is at maximum in a finite system. 
Fig. 2d gives the maximums of S2 around 5.5 MeV for different isotopes, respectively. Again, the critical behavior 
occurs in the same temperature as other observables. 

In conclusion, the critical behaviors are explored for Xe isotopes in the lattice gas model, namely the minimum of 
power-law parameter r of mass distribution, the rise and fall of mean multiplicity of IMF, information entropy and 
Campi's second moment. In a narrow region of critical point, the features of the above quantities show no dependence 
on the isospin of the disassembling system. It reflects that a universal law exists for the same element in the critical 
point. On the contrary, these quantities have isospin dependence at the same temperature outside the critical region. 
Noting that the information entropy is introduced into such an analysis for the first time, and it seems to be useful for 
the searching of critical phenomena in nuclear physics. It will be interesting and meaningful to have some experiments 
to compare our conclusion in the near future. 
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Figure Captions 



Fig.l: Mass distribution of 129 Xe at T=4, 5, 6 and 7 MeV. The lines are the power-law fit. 

Fig. 2: Critical observales: the r parameter from the power law fit to mass distribution (a), the average multiplicity 
of intermediate mass fragments (b), the information entropy (c) and the Campi's second moment (d) as functions of 
temperature and isospin. 
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